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1 Introduction 



rS I The existence, uniqueness and full-rank representation of the Moore-Penrose 

c^ ■ inverse over the quaternion skew field EI have been studied in particular in 



At the same time the problem of determinantal representation of 
the Moore-Penrose inverse of a quaternion matrix is still unsolved. Moreover, 
by means of any introduced determinants currently it is not possible to receive 
a determinantal representation of an inverse matrix over H. The problem 
consists in the fact that there is no determinant functional over a skew field 
which would keep all the properties inherent in the complex case, as is proved 
in [1]. So the determinants of Dieudonne [5] and Study p^ assume values 
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not in a skew field, but in the field which is its center. These determinants 
can not be expanded by cofactors along an arbitrary row or column of a 
matrix. The Moore determinant is introduced p3] in terms of permutations 
only on the class of quaternion Hermitian matrices. The other determinant 
in terms of permutations, the determinant of L. Chen |1] dissatisfies a key 
property of a determinant, its singularity for noninvertible matrices. The 
double determinant introduced by L. Chen can not be expanded by cofactors 
along an arbitrary row or column of a matrix as well. 

In this paper determinantal representations of the Moore-Penrose inverse 
are obtained within the framework of the theory of new matrix functionals 
over the quaternion skew field (the column and row determinants) introduced 
in [TT] . In the first point we cite some provisions from the theory of the col- 
umn and row determinants which are necessary for the following. The theory 
of the column and row determinants of a quaternionic matrix is considered 
completely in [11] . In the second point some known facts from the theory of 
eigenvalues of a quaternion matrix and its singular value decomposition are 
considered. The concept of a characteristic polynomial of a Hermitian matrix 
is introduced and its coefficients are investigated. A Moore-Penrose inverse 
over a quaternion skew field and its limit representation are introduced in 
this point as well. In the third point the theorems about determinantal rep- 
resentations of the Moore-Penrose inverse and the projection matrices A^ A, 
A A"*" for an arbitrary matrix A G H™^" are proved. In the fourth point the 
Cramer rules for a least squares solutions of right and left systems of linear 
equations over a quaternion skew field are obtained. 

2 Elements of the theory of the column and 
row determinants. 

Let M (n, H) be the ring oi n x n quaternion matrices. By H™^" denote 
the set of all m x n matrices over the quaternion skew field H and by HJ?^" 
denote its subset of matrices of rank r. Suppose S'„ is the symmetric group 
on the set /„ = {!,..., n}. 

Definition 2.1 The ith row determinant of A = {a-ij) G M(n,E[) is defined 
by 



rCletjA — / ^ [ i-) ^iiki'^ik-i^iki+l ■ ■ ■'^iki+lii ' ' ' ^ik^ikr + l ' ' ' ^ikr+lr^kr 

uGSn 

for all i = 1, . . . ,n. The elements of the permutation a are indices of each 
monomial. The left-ordered cycle notation of the permutation a is written as 
follows, 

Cr = {i ikiikl+l ■ ■ ■ iki+h) {ik2^k2 + l ■ ■ ■ ^2+^2) • ■ ■ {ikjkr + i ■ ■ ■ ikr + ir) ■ 

The index i opens the first cycle from the left and other cycles satisfy the 
following conditions, ik^ < ikz < ■ ■ ■ < ikr 0'''^d ikt < ikt+s for all t = 2, . . . ,r 
and s = 1, . . . ,lf 

Definition 2.2 The jth column determinant of A = {uij) G M(n,E[) is 
defined by 



Cdet^A= J2 (-1)" ''^IkrJkr+lr 

reSn 



■ '^jkr + likr ■ ■ -^jjki+li ■ ■ ■ ^jki+ljki^jkij 



for all j = l,...,n. The right-ordered cycle notation of the permutation 
T E Sn is written as follows, 

'T = ijkr + lr ■ ■ -jkr + ljkr) ■ ■ ■ {jk2+l2 ■ ■ ■ j k2 + lj k2) {jki+h ' • -jkl + ljkj) ■ 

The index j opens the first cycle from the right and other cycles satisfy the 
following conditions, jk2 < jk-j, < ■ ■ ■ < Jk o.'iT'd jkt < jkt+s for all t = 2, . . . ,r 
and s = 1, . . . ,lt. 

Suppose A*-^ denotes the submatrix of A obtained by deleting botli tlie 
ith row and the jth column. Let a ^ be the jth column and aj be the ith 
row of A. Suppose A ,,• (b) denotes the matrix obtained from A by replacing 
its jth column with the column b, and Aj. (b) denotes the matrix obtained 
from A by replacing its zth row with the row b. 

We note some properties of column and row determinants of a quaternion 
matrix A = (ajj), where i e /„, j G J„ and /„ = J„ = {1, . . . , n}. 

Proposition 2.1 fll^ If b e M, then rdetjAj. {b ■ aj.) = b ■ rdetjA for all 
i = 1, . . . ,n. 



Proposition 2.2 /TT]/ If b e B., then cdetj A.j (a^fe) = cdet^ A& for all j = 
l,...,n. 

Proposition 2.3 01] / ///or A G M(n,EI) it/iere ea;zsts t e In such that 
atj = bj + Cj for all j = 1, . . . ,n, then 

rdetj A = rdet, A^. (b) + rdetj A^. (c) , 
cdetj A = cdetj Aj. (b) + cdetj Aj. (c) , 

where b = (6i, . . . , bn), c = (ci, . . . , c„) and for all i = 1, . . . ,n. 

Proposition 2.4 [77] / If for A G M(n,EI) t/iere exisfo t E Jn such that 
o-it = bi + Ci for all i = 1, . . . ,n, then 

rdetj A = rdetj A . < (b) + rdetj A . t (c) , 
cdetj A = cdetj A . t (b) + cdetj A . t (c) , 

where b = (6i, . . . , bn)^ , c = (ci, . . . , c„)-^ and for all j = 1, . . . , n. 

The following lemmas enable us to expand rdetj A by cofactors along the 
ith row and cdetj A along the jth column respectively for all i,j = l,...,n. 

Lemma 2.1 [Tl^ Let Rij he the right ij-th cofactor of A E M(n,HI), that 

n 

is, rdetj A = ^ Ojj ■ Rij for all i = 1, . . . ,n. Then 



^ f -rdetj A^i (a ,) , z ^ j, 
^^^ \ rdetfcA", t = j, 



where A" (a. j) is obtained from A by replacing the jth column with the ith 
column, and then by deleting both the ith row and column, k = min {/„ \ {i}}. 

Lemma 2.2 /Ti] / Let Lij be the left ij-th cofactor o/ A G M (?7,,HI), that is, 

n 

cdetj A = ^ Ljj ■ ttjj for all j = 1, . . . ,n. Then 



i=l 



^ J -cdetj A[-'(aj.), ij^j, 
'^ ^ cdetfcA-'^ i = j, 



where Af (aj.) is obtained from A by replacing the ith row with the jth row, 
and then by deleting both the jth row and column, k = min{J„ \ {j}}. 



We recall some well-known definitions. The conjugate of a quaternion a = 
Oq + ttii + a2J + 03/? G EI is defined by a = ao — aii — a2J — a^k. The Hermitian 



adjoint matrix of A 



G H'^^'^ is called the matrix A* = (a*,) 



a*j = aji for all i 



m. 



The matrix A 



'^3 ^ mxn 

G H"^"* 



1, . . . , n and j = 1, 
is Hermitian if A* = A. 

A following theorem has a key value in the theory of the column and row 
determinants. 



Theorem 2.1 fH] If A 

cdetiA = 



rdet„A 



j) G M{n,M.) is Hermitian, then rdetiA 
cdet„A G M. 



Taking into account Theorem 12.11 we define the determinant of a Hermitian 
matrix by putting det A := rdetjA = cdetjA for all i = l,...,n. This 
determinant of a Hermitian matrix coincides with the Moore determinant. 
The properties of the determinant of a Hermitian matrix are considered in 
[TT] by means of the column and row determinants. Among them we note 
the following. 

Theorem 2.2 JTl]] If the ith row of a Hermitian matrix A G M(n,EI) is 
replaced with a left linear combination of its other rows, i.e. a^. = Ciaj^. + 
• • • + CkSii^,, where q G H for all I = 1, . . . ,k and {i, ii} C /„, then 



cdetj Aj. (ci • aj, + 



Ck ■ aj. 



rdeti Aj. (ci ■ aj, H \- Ck ■ a^ 



0. 



Theorem 2.3 /ii/ // the jth column of a Hermitian matrix A G M {n,M.) 
is replaced with a right linear combination of its other columns, i.e. a^ = 
a jjCi + ■ ■ ■ + SLj^Ck, where q G H for all I = 1, . . . ,k and {j,ji} C Jn, then 






cdetjA.i (a.j^ ■ ci H h a.^^ ■ c^) = rdet^A. 



ci H ha.j^ ■ Ck) 



0. 



The following theorem about determinantal representation of an inverse ma- 
trix of Hermitian follows immediately from these properties. 

Theorem 2.4 [11] If a Hermitian matrix A G M (n, H) is such that det A 7^ 
0, then there exist a unique right inverse matrix {RA)~^ and a unique left 
inverse matrix {LA)~^, and {RA)~ = (LA)~ =: A~^. They possess the 
following determinantal representations : 

I Rll R2I ■ ■ ■ Rnl\ 



(RA) 



-1 



1 
det A 



R 

R^2 R 



21 
22 



R. 



n2 



\Rln 



R 



2n 



Rr, 



(LA) 



-1 



det A 



L12 

\Lln 



L2I 
L22 

L2n 



Lnl 
Ln2 



\ 



^J 



Here Rij, Lij are right and left ij-th cof actors of A respectively for all i,j = 
l,...,n. 

To obtain determinantal representation of an arbitrary inverse matrix A~^, 
we consider the right AA* and left A*A corresponding Hermitian matrix. 

Theorem 2.5 11 /]/ // an arbitrary column of A E M"^^" is a right linear 
combination of its other columns, or an arbitrary row of A* is a left linear 
combination of its others, then det A*A = 0. 

Since the principal submatrices of a Hermitian matrix are Hermitian, the 
principal minor may be defined as the determinant of its principal submatrix 
by analogy to the commutative case. We introduce the rank by principal 
minors that is the maximal order of a nonzero principal minor of a Hermitian 
matrix. The following theorem determines a relationship between it and the 
rank of a matrix defining as ceiling amount of right-linearly independent 
columns or left-linearly independent rows which form basis. 

Theorem 2.6 |7i]/ A rank by principal minors of A* A is equal to its rank 
and a rank of A e H™-^". 

Theorem 2.7 pJj If A E jjmxn.^ then an arbitrary column of A is a right 
linear combination of its basis columns or an arbitrary row of A is a left 
linear combination of its basis rows. 

The criterion of singularity of a Hermitian matrix is obtained. 

Theorem 2.8 [TTJ The right-linearly independence of columns of A G M™'^" 
or the left-linearly independence of rows of A* is the necessary and sufficient 
condition for det A* A ^ 0. 

Theorem 2.9 JHy // A G M {n, M), then det AA* = det A*A. 

A concept of the double determinant is introduced by this theorem. This 
concept was initially introduced by L. Chen in [4j. 

Definition 2.3 The determinant of the corresponding Hermitian matrix of 
A G M(n,H) is called its double determinant, i.e. ddetA := det (A* A) = 
det (A A*). 
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The relationship between the double determinant and the noncommu- 
tative determinants of E. Moore, E. Study and J. Diedonne is obtained, 
ddetA = Mdet (A* A) = SdetA = Ddet^A. But unlike those, the double 
determinant can be expanded along an arbitrary row or column by means of 
the column and row determinants. 

Definition 2.4 Suppose A G M(n, H). We have a column expansion of 
ddetA along the jth column, ddetA = cdetj (A* A) = X]^«i ' '^ij' '^'^^ ^ ''"'^'^ 

i 

expansion of it along the ith row, ddetA = rdetj (AA*) = ^ aij-M.ij for all 

j 
i,j = l,...,n. Then by definition of the left double ijth cofactor we put hij 

and by definition of the right double ijth cofactor we put Mjj. 



Theorem 2.10 ni] / The necessary and sufficient condition of invertibility 
of Ae M(n,e) is ddetA ^ 0. Then there exists A"^ = (LA)~^ = (i?A)~\ 
where 



(LA)-' = (A*A)^'A* 



1 



/Ln L21 

IL112 1^22 



ddetA 



/Mil 



L^ 



{RA) 



-1 



A* (AA' 



,-1 



1 



ddetA 



V 



12 



In 



'2ra 



21 
22 



2n 



L„i\ 

IL'n2 



.A 



■n2 



(2) 
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and'Lij = cdetj(A*A).j (a* 



rdeti(AA*)j. (a*) for all i,j 



l,...,n. 



This theorem introduces the determinantal representations of an inverse 
matrix by the left ([1]) and right (|2]) double cofactors. 

3 The singular value decomposition and the 
Moore- Penrose inverse of a quaternion ma- 
trix 



Due to the noncommutativity of quaternions, there are two types of eigen- 
values. 



Definition 3.1 Let A G M(n,HI). A quaternion A is said to be a right 
eigenvalue of A i/ A • x = x ■ A for some nonzero quaternion column-vector 
X. Similarly X is a left eigenvalue z/ A ■ x = A ■ x. 

The theory on the left eigenvalues of quaternion matrices has been inves- 
tigated in particular in P [ISl IB]- The theory on the right eigenvalues of 
quaternion matrices is more developed. In particular we note [21El[H]. From 
this theory we cite the following propositions. 

Proposition 3.1 I21^ Let A G M{n,M.) is Hermitian. Then A has exactly 
n real right eigenvalues. 

Definition 3.2 Suppose U G M (n, H). //U*U = UU* = I, then the matrix 
U is called unitary. 

Proposition 3.2 JH]/ Let A G M(?7,,EI) he given. Then, A is Hermitian 
if and only if there are a unitary matrix U G M (n, H) and a real diagonal 
matrix D = diag (Ai, A2, . . . , A„) such that A = UDU*, where Ai, ..., A„ are 
right eigenvalues of A. 

Suppose A G M (n, H) is Hermitian and A G M is a right eigenvalue, then 
A-x = x-A = A-x. This means that all right eigenvalues of a Hermitian 
matrix are its left eigenvalues as well. For real left eigenvalues, A G M, the 
matrix AI — A is Hermitian. 

Definition 3.3 // t G R, then for a Hermitian matrix A the polynomial 
Pa (t) = det (tl — A) is said to be the characteristic polynomial of A. 

The roots of the characteristic polynomial of a Hermitian matrix are its real 
left eigenvalues, which are its right eigenvalues as well. We shall investigate 
coefficients of the characteristic polynomial as in the commutative case (see, 
e.g. [I3j). At first we prove the auxiliary lemma. 

Lemma 3.1 Let A G M (n,H) be Hermitian and the columns ii,. . . ,ik of 
A coincide with the unit vectors ej^ , . . . , ej^ . Then det A equals a principal 
minor obtained from A by deleting the rows and columns ii, . . . ,ik. 

Proof . We claim that if A G M {n, H) is Hermitian and the columns ii, . . . ,ik 
of A coincide with the unit column vectors e^^ , . . . , ej^ respectively, then the 
rows ii, . . . ,ik coincide with the unit row vectors e^^, . . . , ej^ as well. Using 
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Lemma 12.11 we expand det A along the i ith column, where a^^ ^ = for all 
k ^ ii and aj^j^ = 1. Then we obtain 

det A = cdetjj A = 

= -cdeti^A^^^ (ai.) ■ + ■ ■ ■ + cdetiA*i*i ■ 1 + cdet^.A,';" (a„.) ■ = 

= cdetiA*i*i. 

Since the submatrix A*^*^ is obtained from A by deleting both the Zi-th rows 
and columns, by Theorem 12. II it follows that cdetiA'^*^ = det A*^*!. Now we 
calculate this principal minor expanding along the i2-th column. Similarly 
to above we have that det A is equal to a principal minor obtained from A 
by deleting both the zith and i2th rows and columns. Continuing this line of 
reasoning we complete the proof of the lemma. ■ 

Taking into account Lemma 13.11 we can prove the following theorem by 
analogy to the commutative case (see, e.g. [T3]). 

Theorem 3.1 If A e M(n,e) is Hermitian, then pA{t) = T - diP~^ + 
(i2^"~^ — . . . + (—1)" dn, where d^ is the sum of principle minors of A of order 
r, 1 < r < n, and dn = det A. 

For any quaternion matrix A G M{n,M), the eigenvalues of A* A are all 
nonnegative real numbers [17j . 

Definition 3.4 Let A G H'"^". The nonnegative square roots of the n eigen- 
values of A* A are called the singular values of A. 

A key value for a determinantal representation of the Moore- Penrose inverse 
over the quaternion skew field is the following singular value decomposition 
(SVD) theorem. 



Theorem 3.2 [Til (H (SVD) Let A G H™^". Then there exist unitary 



quaternion matrices Ui G M™^™ and U2 G H"^" such that 



U1AU2 



D^ 




G e"^^", (3) 



where D^ = diag (ai, a2, ■ ■ ■ , cXr) , o"i > a2 > ■ ■ ■ > cr.r > 0, and cri, a2, . . . , (Tj. 
are the all nonzero singular values of A. 



As unitary matrices are invertible, the equality (E]) can be written as follows 

A = VEW*, (4) 

where V G H'"^"* and W G H"^" are unitary matrices, and the matrix S = 

{(Tij) G H™^" is such that an > 0"22 > ■ ■ ■ > Urr > CTr+lr+l = ■ ■ ■ = aqg = 0, 

q = min{n, m}. 

We get the following lemmas, which have the analogues in the complex 
case [3j. 

Lemma 3.2 Suppose A G H™"^" has the singular value decomposition, A = 
VSW*. Let A+ = W ■ S+ ■ V*, where S+ G H"^*" is obtained from E by 
transposition and replacing positive entries of Xl by reciprocal. Then for A"*" 
the following conditions are true 



1) (AA+)* = AA+; 

2) (A+A)* = A+A; 

3) AA+A = A; 

4) A+AA+ = A+. 



(5) 



Proof. We obviously have {^^) = S and ( (S+) j = S"*" for S from 
the SVD by (g]) and S+. Then it follows that 

(AA+)* = (VSW*WS+V*)* = (VSIS+V*)* = Tv (S+)^ S^V*)* = 
= ('V(S+)^ W*W5]^V*y = V5]W*WS+V* = AA+. 

The proof of 1) is completed. By analogy we can prove 2). 

Now we prove the condition 3). Note that SS+ = "■ ^ e H™^™. 

This implies SS+S = S, then AA+A = VSW* ■ WE+V* ■ VSW* = 
V ■ SS+S • W* = V ■ S ■ W* = A. 

By analogy to 3) can be prove the condition 4).B 

Lemma 3.3 There exists a unique matrix A+ that satisfies conditions l)-4) 
in ^. 
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Proof. Suppose that both matrices B G H"'^'^ and C G H"-^™ satisfy condi- 
tions l)-4) in Lemma [3.21 Then we have 

B = BAB = B (AB)* = BB*A* = BB* (ACA)* = BB*A*C*A* = 

= B (AB)* (AC)* = BABAC = BAG = BACAC = (BA)* (CA)* C = 
= A*B*A*C*C = (ABA)* C*C = A*C*C = (CA)* C = CAC = C. 



Definition 3.5 Let A G H'"^". The matrix A"*" is called the Moore-Penrose 
inverse if it satisfies all conditions in ^. 

By analogy to the complex case ^ we have the theorem about the limit 
representation of the Moore-Penrose inverse. 

Theorem 3.3 //AG H™^"' and A"*" is its Moore-Penrose inverse, then 
A+ = limA* (AA* + al)"^ = lim (A* A + al)~^ A*, where a G IR+. 

a— s>0 a— s>0 

Proof. Suppose A = VSW*, then A* = WS*V* and A+ = WS+V*. 
Since V is unitary, then V* = V^^. We have 

A* (AA* + aiy^ = WSV* • (VS ■ W*W • SV* + al)~^ = 
= WSV* • (V (SS* + al) V*)"^ = WS (SS* + al)"^ V*. 

Consider the matrix 



/ 



S (SS* + al) 











\ 







A„- 



V 







/ 



It is obviously that limS (SS* + al) ^ = S+. This implies lim A* (AA* + al) ^ 

cn-O a— >-0 

limWS (SS* + aiy^ V* = A+. 

By analogy we can prove that A+ = lim (A*A + al)^ A*.B 

a— >0 

Corollary 3.1 //AG H™^"', then the following statements are true, 
i) //rank A = n, then A+ = (A*A)"^ A* . 

a) //rank A = m, then A+ = A* (AA*)"^ . 

Hi) //rank A = n = m, then A+ = A^-*^ . 
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4 Determinantal representation of the Moore- 
Penrose inverse. 

Lemma 4.1 If A e H™^", then rank (A* A) . (a*,) < r. 

Proof. Let's lead elementary transformations of the matrix (A* A) ^ (a*) 
right-multiplying it by elementary unimodular matrices Pj^ (—0,^)5 ^ ¥" 3- 
The matrix Pj^ ("ti ;,) has —ajk in the (z, A;) entry, 1 in all diagonal entries, 
and in others. It is the matrix of an elementary transformation. Right- 
multiplying a matrix by Pj^ (—0,^)5 where k 7^ j, means adding to A;-th 
column its z-th column right-multiplying on — a^fc. Then we get 



(A*A).,(a*^.)-nP^^'(-«^-^) 



/ Y. ^ifc^fci • • • a\j 



Yl O'lk'^kn \ 



\ Y. Kk'^kl ■ ■ ■ <j • • • Y Kk'^kr 
\ k^j k^j 






The obtained matrix has the following factorization. 



/ Y C'lk'^kl ■ ■ ■ alj ... Y 0*lk^kn \ 
k^j k^j 



\ Y K.k'^ki 
\ k^j 



'^■21 0^22 



Denote by A : = 







\ CLml 



* \ 



^2m 



\ ^nl ^n2 • • • ^nm / 

f ail ... 



n ■ 
i—th 

( an 




Y '^nk^kn 
ky^j 





1 





i—th 



inl \ 







J - th. 







'th 



dnl \ 



j — th. The matrix A is ob- 



(^mn J 



tained from A by replacing all entries of the jth row and of the ith column 
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with zeroes except that the {j,i) entry equals 1. Elementary transforma- 
tions of a matrix do not change its rank and the rank of a matrix product 
does not exceed a rank of each factor. It follows that rank (A* A) ■ (a*) < 

min < rank A* , rank A >. It is obvious that rank A > rank A = rank A*. 
Taking into account Theorem 12.61 we obtain rank A* A = rank A. This com- 
pletes the proof. ■ 

The following lemma is proved in the same way. 

Lemma 4.2 If A e H™^", then rank(AA*) . (a*,) < r. 

We shall use the following notations. Let a := {ai, . . . , a^} C {1, . . . , m} and 
/3 := {/3i, . . . , f^k} C {1, . . . , n} be subsets of the order 1 < A; < min {m, n}. 
By A'a denote the submatrix of A determined by the rows indexed by a and 
the columns indexed by /?. Then A" denotes the principal submatrix deter- 
mined by the rows and columns indexed by a. If A G M (n, H) is Hermitian, 
then by | A^l denote the corresponding principal minor of det A. For 1 < A; < 
n, denote by Lk^n := { a : a = (ai, . . . , a^) , 1 < cti < ■ ■ ■ < ctfc < n}, the col- 
lection of strictly increasing sequences of k integers chosen from {1, . . . , n}. 
For fixed i E a and j G f3, let Ir m{i} '■= {a : a E Lrm,i ^ «}, Jr n{j} '■ = 
{/3:/3gL,,„,jG/3}. 

Lemma 4.3 //AG H"^^" and t G R, then 

cdet, (tl + A*A),^ (a*,.) = cf^'^t^^^ + c^^'k^-^ + ■■■ + c^:^\ (6) 

where c^^^ = cdet^ (A*A) ,. (a*^.) and 4'^'^ = ^ cdet^ ((A*A) . (a*^.)) ^ 

/3eJfc,„{j} 
for all k = 1, . . . ,n — 1, i = 1, . . . ,n, and j = 1, . . . ,m. 

Proof. Denote by b j the i-th column of the Hermitian matrix A*A =: 
i^v)nxn- Consider the Hermitian matrix (tl + A* A) • (b. j) G H"^". It differs 
from (tl + A* A) by an entry bu. Taking into account Theorem 13. II we obtain 



det (tl + A*A) . (b. ,) = dir-i + d2t"~' + --- + dn, (7) 
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where dk = Yl ^^^ (A* A) ^ is the sum of all principal minors of order 



/3eJfc,„{i} 

k that contain the z-th column for all A; = 1, . . . , n — 1 and dn = det (A* A) 
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Consequently we have b. 



I 






Y^a*iaii, where a*; is the Ith 

I 



J 



column- vector of A* for all/ = 1, . . . , m. Taking into account Theorem 12 ■![ 
Lemma 12.11 and Proposition 12.21 we obtain on the one hand 

det (tl + A* A) . (b. i) = cdet^ (tl + A* A) . (b, i) = 
= Yl cdet, (tl + A* A)', {a*,au) = ^ cdet^ (tl + A*A) , (a*J ■ a^ (§) 



On the other hand having changed the order of summation, we get for all 

k = 1, . . . ,n — 1 

4= E det (A* A) ^= E cdet,(A*A)^ = 
/3eA-,n{i} /3eJfe.„{i} 

E Ecdet, ((A*A). Ja*,aH))? = E E cdet, ((A*A)^ (a:,)) ^ ■ a^. 

/3eJfe,„{i} I I l3eJk.n{i} 

(9) 
By substituting ([8]) and ([9]) in ([7]), and equating factors at an when / = j, 
we obtain the equality (|6]). ■ 

By analogy the following lemma can be proved. 

Lemma 4.4 If A e H™^" and t e R, then 

rdet,(tl + AA*),-. (a*) = rf'k''^^ + rfh""-^ + ■■■ + r^^^l 



where r^n^^ = rdetj(AA*)j, (a* ) and rf^ = E rdet^- ((AA*)^-. (a* )) 

aelr,m{j} 

for all k = 1, . . . ,n — 1, i = 1, . . . ,n, and j = 1, . . . ,m. 



Theorem 4.1 //AG H^^", then the Moore-Penrose inverse A+ = (0^*^) G 
^nxm pQgg^gg ifi^ following determinantal representations : 



I l3eJr,n.{i} 

4 = 



E cdet,((A*A).,(a^)) 



JJJ /9 



E 

/3GJr, n 



'A*A) 



(10) 
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or 



E rdet,((AA*),.(aI)) 

_|_ "G/r,m{j} 



*J 



E l(AA*)2| 



(11) 



Proof. At first we prove (TTO]) . Using Theorem IO| we get A+ = lim (al + A* A) ^ A*. 

a— !>0 

The matrix (al + A* A) G M"^" is a full-rank Hermitian matrix. Taking into 
account Theorem 12.41 it has an inverse, which we represent as a left inverse 
matrix 



(aI + A*A)"^ 



/ Lii L21 . . . Lni \ 
L12 L22 ■ ■ ■ Ln2 



det {al + A* A) 



\ Lin L 



2n 



^nn / 



where Lj,- is a left zjth cofactor of a matrix al + A*A. Then we have 



-1 



(aI + A*A)"'A 



/ ra n 

I Y^ -^fcl«fcl Yj LklO*k2 



1 



det(aI+A*A) 



fc=l 



fc=l 



v 



Z] -^fc2«fci Z] -^fc2afc2 
fc=l fc=l 



Zj LknO*kl zJ LknO*k2 
k=l k=l 



* \ 



Zl Lkiakm 
fc=i 

n 

Z] Lk2alm 
k=l 



Z^ Lkn^km t 
k=l / 



Using the definition of a left cofactor, we obtain 



A^ 



lim 



/ cdeti(Q!l+A*A) i(a*j] 
det(aI+A*A) 



V 



cdct„(Ql+A*A)^(a*^j 
det(Q:I+A*A) 



cdeti(aI+A*A) ]^(a*„J \ 

dct(aI+A*A) ' 

cdctn(aI+A*A),^(a*^) . 

det(aI+A*A) / 



where dk = Z^ 



'A*A) 



(12) 



By Theorem Owe have det (al + A* A) = a" + rfia""^ + ^20""^ H h 4 



is a sum of principal minors of A* A of order k 



for all k = 1, . . . ,n — 1 and dn = det A* A. Since rank A* A = rank A = r 
and dn = dn-i = ■ ■ ■ = dr+i = 0, it follows that det (al + A* A) = a" + 
dia"-i + ^2"""^ + ■ ■ ■ + 4a"-^ Using (0) we get cdet^ (al + A*A) . (a*,.) = 
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Ci a"' ^ + C2 a"' ^ + ■ ■ ■ + Ok'' for all i = 1, . . . , n and j = 1, . . . ,m, where 
Cfc^^ = E cdeti ((A*A).i (a*,.)) ^ for all A; = 1, ... ,n - 1 and cH^^ = 

l3eJk,n{i} 

cdet,(A*A) Ja*,.). 

Now we prove that c^*"' = 0, when k > r + 1 for all z = l,...,n, 
and j = 1, . . . ,m. By Lemma 14.11 rank (A* A) ^ (a*) < r, then the matrix 
(A* A) ■ (a* ) has no more r right-linearly independent columns. 

Consider ((A* A) j (a*,)) ^, when /? G Jfc.nIO- It is a principal submatrix 
of (A* A) ^ (a*) of order k > r + 1. Deleting both its z-th row and column, 
we obtain a principal submatrix of order fc — 1 of A* A. We denote it by M. 
The following cases are possible. 

Let k = r + 1 and det M 7^ 0. In this case all columns of M are 
right-linearly independent. The addition of all of them on one coordinate to 
columns of ((A* A) ■ (a*)) ^ keeps their right-linear independence. Hence, 
they are basis in a matrix ((A* A) ■ (a*)) ^, and by Theorem 12.71 the i-th. 
column is the right linear combination of its basis columns. From this by 
Theorem 12. 5t we get cdeti ((A* A) ■ (a*j)) ^ = 0, when /3 G Jk,n{'>'} and 
k>r + l. 

If A; = r + 1 and det M = 0, then p, {p < k), columns are basis in 
M and in ((A*A) . (a*,)) ^. Then by Theorems O and ES] we obtain 

cdeti ((A*A) . (a*,.)) ^ = as well. 

If A; > r + 1, then from Theorems 12.61 and 12.81 it follows that det M = 
and p, {p < k — 1), columns are basis in the both matrices M and 
((A* A) i(a*j)) ^. ThenbyTheoremsOandESlweobtainthatcdeti ((A*A) , (a*)) ^ 
0. ' 

Thus in all cases we have cdetj ((A* A) ■ (a*^)) ^ = 0, when /3 G Jk,n{i} 

and r+1 < k < n. From here if r+1 < k < n, then c]^^' = Yl ^detj ((A* A) . (a*^)) |^ 













/3eJfc,„ 


{i} 




0, 


and Cn 


= cdeti (A*A) . (a*,.) 


= for alH = 1, 


. . . ,n and j 


=1,.... 


m. 




Hence, 


cdetj (al + A*A) . 


(a* 


^.)=cf^)a"-i + c 


2 « H 


+ cl'^')«" 


-^for 


all 


.^ = 1,.. 


, . ,n and j = 1, . . . 


, m. 


By substituting 


these values 


in the matrix 
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from (1121) . we obtain 



/ 



.(ll)^n~l 



(11). 



A^ 



lim 



a"~i + ---+4 'a""'' 



.("l)^n-l_^..._^^(nl)^n-r 



(l7n) T)— 1 



^(l"i)^n-r 






/ .(11) 



V 



dr 

J"l) 



(im) \ 



dr 

(nm) 



J 



Here c 



(ii) 



'A*A) 



E cdet, ((A*A) , (a*,.)) ^ and d. = E 

;SGJr, n{i} /3GJr, n 

Thus, we have obtained the determinantal representation of A"*" ( ITOj) . 

Similarly one can prove the determinantal representation of A+ f fTTl) . ■ 

Remark 4.1 /n / fTOj) the index i in cdetj ((A* A) ^ (a*)) « designates ith 

column o/((A*A) ^ (a*)), but in the submatrix ((A* A) ^ {^*j)) r ^^^ entries 
of a*- may 6e placed in a column with the another index. In / flT]) we have 
equivalently. 

Remark 4.2 //rank A = n, then by Corollary\3l^A~^ = (A*A)"^ A*. Con- 
sidering (A* A) as a left inverse, we get the following representation of A~^: 



'cdeti(A*A).i(a*i^ 



A^ 



ddetA 



cdeti(A*A).i(a*„ 
cdet„(A*A).„(a*^).^ 



(13) 



^cdet„(A*A).„(a*i) . 
If m > n, then by Theorem \4.1\ for A+ we have l[Td\) as well. 



Remark 4.3 //rank A = m, then by Coronary (231 A+ = A* (AA*) \ 
Considering (AA*)^ as a right inverse, we get the following representation 
ofA+: 



A^ 



1 



^rdeti(AA*)i. (a]; 



rdet^(AA*)^. (aj 



ddetA 



rdetm(AA*)^,(a; 



jdeti(AA*)i.(a;.) . 
If m < n, then by Theorem \4.1\ for A+ we also have DJ]) 



(14) 



17 



Corollary 4.1 If A E M.^^"^, where r < min {m, n} or r = m < n, then for 
a projection matrix A+A =: P = {Pij)nxn '^^ have its following determinantal 
representation 

E cdet,((A*A),(d,))? 

/3eJr. n{i} 

P^j = 



E 

l3€Jr,n 



:a*a) 2 



where dj is the j-th column of A* A G H"^" and for all i,j = l,...,n. 

Proof. Representing A^ by (TTOl) and right-multiplying it by A, we obtain 
for an entry pij of A+A =: P = (pij)„x„- 



E^ cdet,((A*A).,(a*^.))^ 

Cip ■ CLpj — 2_^ = \(A*A)'^\ ' ^^-^ ~ 

E Ecdot,((A-A)^(a:^.))^-afe, E cdet.((A*A). ,(d.,)) 

/36Jr, n{i} ^ /Se-Tr, n{i} 



E |(A-A)^| E |(A*A)'^| 

/3SJr, n /3eJr, n 

where d j is the jth column of A* A G H"^" and for all i,j = 1, . . . ,n. ■ 
By analogy can be proved the following corollary. 

Corollary 4.2 If A E ^mxn ^ where r < min {m, n} or r = n < m, then for 
the projection matrix AA"*" =: Q = (%)„xm '^^ have its following determi- 
nantal representation 

E |((AA*)ag..))Sl 
^'^' ^ E l(AA*)-| ' 

where gj, is the jth row of (AA*) G M*"^™ and for all i,j = l,...,m. 

Remark 4.4 By definition of a classical adjoint matrix of A G C"^", Adj [A]- 
A = A • Adj [A] = det A ■ I. Let A G H™^". //rank A = n, the by Corol- 
lary \3.1\ we have A+A = !„. Representing A+ by ( fl^) as A+ = j^^.^,^^ , 

where L = (cdetj ((A*A) . (a*,) ))„^^, we o6iam LA = det (A*A) ■!„. This 
means that the matrix L =: Adj l [A] is the left classical adjoint matrix of 
A G H*"^". 



// rank A = m, then by definition of a right classical adjoint matrix of 
A E M™^" by Corollary \y.l\ and by [T^ we can put 



Adj«[A]:=((rdet,(AA*),,(aD)U„, 

since in this case A ■ Adj r[A] = det(AA*) ■ I. 

//rank A = r < min{m,n}, then an analog of a left classical adjoint 
matrix of A E H™^" by / fi^) can accept 



Adji[A]:= I Yl cdet,((A*A),,a*^.) 

Indeed, since eigenvalues of a projection matrix are only 1 and 0, there exists 
such a unitary matrix U G H"^"' that 

Adji[A].A= E |(A*A)^|-P = 

= E |(A*A)S|-Udiag(l,...,l,0,...,0)U*. 

//rank A = r < min {?tt,,?7,}, then by an analogue of a right classical 
adjoint matrix of A E M™^" by 177]) we can put 



Adin[A]:= I Y. rdet, ((AA*),, (a*) ) 

as there exists such a unitary matrix V E M™^™ that 

A-Adj«[A]= E |(AA*)-|-Q = 

Ct€Jr, m 

= E |(AA*)°|-Vdiag(l,...,l,0,...,0)V*. 



Remark 4.5 If A E C™^" is a matrix with complex entries, then we obtain 
the following analogs of l[TO\) and 07]) . respectively. 



E ((A*A).Ja*,.)) 

^ /SeJr, n{«} 



E 



J^Jr, n 



E (A*A)^ 






;AA*)^,(aI))- 



E l(AA* 
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for all i = 1, . . . ,n and j = 1, . . . ,m. These determinantal representations 
are original in this case as well. It is reflected in [WJ. The analogous result 
is obtained in WW. 



Remark 4.6 The Gram polynomial of a matrix A, 

detilm + ^AA*) = 1 + aiz + ■ ■ ■ + a^z", 

is used for generalizing the Moore-Penrose inverse and the classical Cramer's 
rule to the corresponding undetermined and overdetermined cases over an 
arbitrary field in f^] as well. 

5 Cramer's rule for a least squares solution 
of quaternion system linear equations 

Definition 5.1 Suppose 

A-x = y (15) 

is a right system linear equations over the quaternion skew field M, where 
A G H™^" is the coefficient matrix, y G H'"^^ is a column of constants, and 
X G H"^-'^ is a unknown column. The least square solution of l[T5\) (with the 
least norm) is called the vector x° satisfying 

II x'^ 11= mill < II X II : II A • X — y 11= min || A ■ x — y 

where H" is an n-dimension right quaternion vector space. 

We recall that in the right quaternion vector space H" by definition of the 
inner product of vectors we put (x, y) := y*x = yl ■ Xi + ■ ■ ■ + y^ ■ x„ and 
||x|| := ^/(x, x) is the norm of a vector x G H". By analogy to a complex 
case (see, e.g. [7]) we can prove the following theorem. 

Theorem 5.1 The vector x = A+y is the least square solution of ( 173]) . 

Definition 5.2 Suppose 

x-A = y (16) 

is a left system linear equations over the quaternion skew field H, where 
A G H™^" is the coefficient matrix, y G M^^"' is a row of constants, and 
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:si E M^^"^ is a unknown row. The least square solution of / fl^) (with the least 
norm) is called the vector x*^ satisfying 



mm 



li- A 



mill llx ■ A — yl 



where ™H[ is an m-dimension left quaternion vector space. 

We recall that in the left quaternion vector space '"H by definition of the 
inner product of vectors we can put (x, y) = xy* = Xi -y^ + ■ ■ ■ + Xm ■ Um- 
Then ||x|| := a/(x, x) is the norm of x G ""EI. 

Theorem 5.2 The vector x = y ■ A+ is the least square solution of l[TS\) . 

Theorem 5.3 (i) //rank A = n, then for the least square solution x° = 
(x°, . . . , x1)^ of 173]) ■ we get for all j = 1, . . . ,n 



X, 



cdet,(A*A),^.(f) 



^ ddetA 

where f = A*y. 

(ii) If rank A = k < m < n, then for all j = 1, . . . ,n we have 

E cdet,((A*A).^.(f) 

_ /3eJr,n{j} 



(17) 



X 



P&Jr, ■ 



E (A*A)^ 



(18) 



Proof . i) If rank A = n, then A"*" can be represented by 013p . Denote 
f := A*y. Representing A+y by coordinates we obtain (fT7|) . 
ii) If rank A = k < m < n, then by Theorem 14.11 we represent the matrix 
A+ by (TTOj) . Representing A+y by coordinates we obtain (TT8|) . ■ 



Remark 5.1 In a complex case the following analogs of (11) and pR) are 
obtained respectively in ^12] for all j = 1, . . . ,n, 



^ det (A* A) ^. (f ) ^ 
^ det r A*A) ' ^ 



E 



((A*A) ^. (f) 



E (A*A) 



/3 
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Theorem 5.4 (i) //rank A = m, then for x° = (x°, . . . ,x^) of i[TS\) we 
obtain for all i = 1, . . . ,m 



X 



_ rdet, (AA*) . (z) 
ddetA 



where z = yA*. 

(a) If rank A = k < n < m, then for all i = 1, . . . ,m we have 

E rdet,((AA*)az))2 

Q oe-fr,m{«} 



X, 



E l(AA* 



The proof of this theorem is analogous to that of Theorem [5? 



(19) 



(20) 



Remark 5.2 In a complex case the following analogs of / iTPj) and [20\) re- 
spectively are obtained in 11^ for all i = 1, . . . ,m, 



.WAA.^ . ^ ^ l((AA*)az)) 

_ det (AA*) i. (z) _ aeir,m{i} 
Xi ; ' : r~ ) Xi 



det AA* 



E l(AA* 



6 Example 



Let us consider the left system of linear equations. 

xii + 2x2^ — X3 = z, 

-Xi/C + X2J + X3J = j, 
S^li + ^2 + 2:3/;; = /c, 

Xi + X2A; + x-^i = 1. 



(21) 



z —k j l\ 
The coefficient matrix of the system is the matrix A = | 2z j 1 A; . 

-1 j k ij 

The row of unknown is x = ( Xi X2 X3 ) and the row of constants is 
y=(z j k 1). Then for (12T]) we have x ■ A = y. We obtain 

/-i -2i -l\ 

k -j -j 

-3 1 -k 

\ 1 -k -i) 
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4 2 -i + j - k -4i 

AA* = \ 2 + i-j + k 7 l-2i-j-k\. 

Ai l + 2i+ J + k 4 

Since ddetA = det AA* = rdetiAA* = and 



det (AA*)^^ = rdeti 



4 2-i+ j -k 

2+i-j+k 7 

4:-7-{2-i + j-k)- {2 + i-j + k) = 21^0 



by Lemma [2.61 rank A = 2. We shall represent A+ by fITT]) . 



E |(AA*)S| = det 
'^'^''Hl + 21 + j + k 4 



4 

2+i-j+k 
l-2i- j-k^ 



2-i+j-k 
7 
4 
Ai 4 



+ 



det 



Ai 



42. 



Now we shall calculate Vij = J2 ^^etj ((AA*) ^.(a*)) ^ for alH = 1,4 

ctei2,3{j} 
and j = 1,3. To obtain rn, we consider the matrix 



(AA*)i.(at.: 



Then we have 



-i -2i -1 

2 + i-j + k 7 1 -2i- J - k 

Ai l + 2i + j + k A 



Til = rdeti 



-2i 



+ rdeti 



'^\2 + i-j + k 7 ) ' ^^^"' \Ai A 
= -i-7 - (-2i) ■ {2 + i - j + k) - i ■ A - {-I ■ Ai) = -2 - ?,i - 2j - 2k, 

and so forth. Continuing in the same way, we get 



f-2 -3i-2j -2k 2- 12i + 2j + 2k -3 + 2i + 2j - 2k\ 

l + i + 2j + Qk -2 + 2i-6j -Ak 1 - i - Qj + 2k 

-2-i-6j -k Q-2i + Aj + 2k -1 + 2i + j - 6k 

\ 6 + i + j + 2k -A + 2i-2j -6k 1 - 6i - 2j + k J 

We find the least square solution by means of the matrix method by Theorem 



A+ = — 
42 



x° = y ■ A+ = — (8+ in + 3j-3A;, 12 - 4? - 8j, 11 - 8i + 3j + 3k) . 
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Now we shall find the least square solution of (1211) by means of Cramer's 
ruleby 020]). Wehave z = y ■ A* = (2 + 2i, 3, 2 - 2?) • Since 

/ 2 + 2z 3 2-2z 

(AA*)i.(z)= 2 + 2-j + A; 7 1 - 2i - j - k 

\ 4i l + 2i+j + k 4 

we get 

, / 2 + 22 3\ , /2 + 22 2 - 2A 

rdeti „ . . , „ + rdeti , . , 

o_ ^\2 + i-j + k 7J^ ^\ At 4: J _8 + lli + 3j-3k 

""' " E l(AA*)-| - 42 • 

4 2-i + j - k -Ai 
Since (AA*)2.(z) = I 2 + 2i 3 2-2i 

Ai l + 2i + i + k A 

, A 2-i + j-k\ ( 3 2-2A 

rdet2 +rdcti 

o_ V2 + 2Z 3 I \l + 2i + 3 + k A j _ 12-4.-8. 

^2- Y. l(AA*)g| - 42 • 

QS-f2, 3 

Since 

4 2- i+ i -k -Ai 

(AA*)3.(z)= 12 + 2-j + A: 7 1 - 2^ - j - A: | , 

2 + 2i 3 2 - 2i 



, 4 -Ai \ (7 1 - 2i - j - k 

rdet2 +rdeti 

^0 _ lv2 + 2z 2-2z; 1,3 2-22 y _ n-8^+3,+3fc 

•^3 •£ l(AA*)g| 42 

a!el2^3 
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